Coupled diabatic potential energy surfaces for studying the nonadiabatic dynamics at conical intersections in angular resolved photodetachment simulations of OHF--->OHF+e-.

An energy-based method is proposed for the diabatization of the OH(2Pi)+F(2P)-->O(3P)+HF(1Sigma+) reaction. It is demonstrated that the diabatic representation obtained is regularized, i.e., the residual derivative couplings do not present singularities at the conical intersections appearing along the reaction path. This method only requires the knowledge of the 1,2 3A" and 1 3A' eigenvalues and does not require any adjustable parameter. Thus, many convergence problems arising in other derivative-based diabatization methods are avoided, and the description of the configuration space along the reaction path is enormously simplified. Three-dimensional coupled diabatic energy surfaces are obtained by an interpolation procedure using approximately 4000 accurate ab initio points. The angular resolved photodetachment cross sections are obtained in the diabatic and adiabatic representations using a wave packet method. An excellent agreement is obtained with recent experimental data [D. M. Neumark, Phys. Chem. Chem. Phys. 7, 433 (2005)] for high electron kinetic energies where only the triplet electronic states contribute.